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a ' Abstract 

We evaluate the one-loop prefactor in the false vacuum decay rate in a the- 
ory of a self interacting scalar field in 3 + 1 dimensions. We use a numerical 
method, established some time ago, which is based on a well-known theo- 
rem on functional determinants. The proper handling of zero modes and of 
renormalization is discussed. The numerical results in particular show that 
quantum corrections become smaller away from the thin-wall case. In the 
thin-wall limit the numerical results are found to join into those obtained by 
a gradient expansion. 
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1 Introduction 



First-order phase transitions play an important role in various phenomena 
from solid state physics to cosmology. The basic theoretical concepts of 
these transitions have been developed long ago [U El El 111 El Ej ■ The phase 
transition proceeds via formation of stable phase (or true vacuum) bubbles 
within a metastable (or false vacuum) environment, and subsequent growth 
of these bubbles. Two mechanisms of the first order phase transitions are 
known: quantum tunnelling and thermal activation. In both cases the decay 
rate of a metastable state is given by the formula 



For tunnelling in a (3 + 1) dimensional theory the quantity B in the exponent 
is given by the classical 4d Euclidean action evaluated on a bounce, a finite 
action Euclidean solution of classical equations of motion, asymptotically ap- 
proaching the false vacuum. For thermal activation at nonzero temperature 
T the exponent is given by —B = S/T, where S is the energy of a critical 
bubble (sphaleron), which is a static solution "sitting" on a top of a barrier 
separating two vacua. The bounce as well as the sphaleron are unstable so- 
lutions with just one negative mode. Bubbles smaller than critical collapse, 
and the ones bigger than critical expand and lead to the transition to a new 
phase. These static solutions and Euclidean solutions are related, namely the 
sphaleron in {d + 1) dimensions can be viewed as a bounce in d dimensions. 

The leading order estimate for the transition rate is easy to obtain, it 
just requires solving - in general numerically - an ordinary, though nonlinear 
differential equation. Analytic estimates can be obtained in the so-called 
thin- wall approximation. 

The pre-exponential factor A in Eq. (jl.lj) is calculated taking into account 
quadratic fluctuations about the classical solution and is given as a ratio of 
the functional determinants. In general it is a very difficult task to calculate 
analytically the determinants, while the background solution itself is not 
known in a closed form. It has taken two decades until the first (numerical) 
computations of the quantum corrections to the leading order semiclassical 
transition rates have appeared [711110 UHl HH Of course nowadays the 
CPU time requirements for such computations are, even for more involved 
systems, of the order of seconds. On the other hand the requirements of a 
precise renormalization, which compares exactly to the one of perturbative 
quantum field theory, and of the inclusion and careful treatment of high 
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partial waves, have of course remained the same. The method used here has 
been developed and tested for various systems and has become a standard 
procedure. It is well suited for computations of coupled channel problems as 
well 1131 . 

While the special technique used here applies only to the computation 
of functional determinants, the general approach can be used as well for 
computing zero point energies ^| via Euclidean Green functions. Of 
course functional determinants can be computed likewise using Euclidean 
Green functions [121 CHI- Various other techniques for computing the exact 
quantum corrections have been developed in the past decade. In Refs. [TT| [T7] 
the heat kernel is computed using a discretisation of spectra, in Ref. [TS] 
Minkowskian instead of Euclidean Green functions are used, and in Ref. ^H] 
the zero point energy is computed via the ( function. 

The effective action may be computed approximatively by using gradient 
expansions. There is an ample literature on this subject. We just quote 
Refs. |2ni HH 1221 12B| for expansions using advanced heat kernel techniques, 
and Ref. P3] for expansions based on Feynman graphs. 

The leading quantum corrections, being essentially a one loop effect, can 
be viewed as a "summary" of the particle creation during the phase transition 
|2Sj. The question about the quantum corrections is very important one, 
while there are cases when particle creation is so strong that it drastically 
modifies the original classical tunnelling solution [221 123 • 

The aim of the present paper is to calculate the pre-factor A for tunnelling 
transitions in a theory of one self- interacting scalar field theory in (3+1) 
dimensions. 

The rest of this paper is organized as follows: In the next Section we will 
describe our strategy for calculation of one loop effective action. In Section |21 
we formulate our model, specify the form of the potential, write the equation 
of motion for the bounce and present our numerical results for classical action 
S[ip]. In Section|3]we describe the calculation of the fluctuation determinant, 
Eq. ()2.3p . There we also discuss regularization and renormalization. Our 
numerical results are presented and discussed in Sectional We end with some 
general remarks and conclusions in Section [HI Formulas describing the thin- 
wall approximation and gradient expansion are collected in the Appendixes 
A and B respectively. 
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2 General strategy 



We will consider phase transitions in a theory of one self interacting scalar 
field (/? in 3 + 1 dimensions. Corresponding Euclidean action is 



where the field potential U {(f) is assumed to have two non-degenerate minima 
ip = if- and (fi = > (compare Fig. ^ and it will be given explicitly in 
the next section. For convenience we have fixed the value of ip in the unstable 
vacuum as V9_ = 0. 

Any state built on the local minimum (p_ is metastable. It can tun- 
nel locally towards the phase. The tunnelling rate per unit volume per 
unit time, 7 = T/VT, is supposed to be dominated by the classical action 
Sci of a field configuration, the bounce (pb{x), which looks like a bubble of 
the v'+'phase within the phase. In particular it can be shown [2H] that 
the bounce configuration ipb{x) which minimizes the action is spherically 
symmetric in four- dimensional Euclidean space. In the tree level approxi- 
mation the decay rate is determined essentially by the tunnelling coefficient, 
7 oc exp{-Sci[ipb{x)]Y. 

The tree level tunnelling rate receives corrections in higher orders of 
the semiclassical approximation. In quantum field theory the fluctuations 
around the bounce contribute in the next-to-leading order approximation a 
pre-exponential factor to the decay rate. The rate per volume and time is 
known to take the form [5] 



The prime in the determinant implies omitting of the four translation zero 
modes. With the second equation we have introduced the fluctuation oper- 
ator in the background of the bounce 




(2.1) 




(2.2) 




(2.3) 



M 



{d/drf - A + U"{^) 



(2.4) 



■^For a more concise statement see Sectional 
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and its counterpart Ai^^^ in the unstable vacuum. 

The counterterm action Set is necessary in order to absorb the divergences 
of the one-loop effective action 

stiLM = l^^\'^M- (2-5) 

In order to evaluate the one loop effective action we decompose fluctua- 
tions about the bounce ipb into 0(4) spherical harmonics, calculate the ratio 
of determinants Ji of partial wave fluctuation operators and obtain lnT> as 
J2idilnJi, where di is the 0(4) degeneracy di = {I + 1)^ (see e.g. I^H])- In- 
calculating InV we exclude the divergent perturbative contributions of first 
and second order in the external field of the bounce (fb- The regularized 
values of these contributions are then added analytically. All divergences of 
InV appear in the standard tadpole and fish diagrams. We will not specify 
Set explicitly, we will equivalently omit the divergent parts of lnP[(/)] using 
the MS convention. 



3 The Tree-Level Action 

In this section we specify our model, discuss the bounce solution and proper- 
ties of corresponding classical action. We parameterize the 99^-potential with 
two minima as 

f/(^) = ^mV-^^' + ^A^S (3.1) 

and choose the same dimensionless variables as in Ref. (SUlEl: = X'^/m 
for /i = 0, 1, 2, 3, and ip = The classical action then takes the form 

Seiiif) = f3Sei{if) , (3.2) 
where rescaled classical action Sci{f) is 

Scii^) = J rf'X Q (V<l>) V t/(<I>)) , (3.3) 

with 

[/($) = - + , (3.4) 



and a and f3 

8 = 

At] 



2 

m 



/3 = — , a = A/3. (3.5) 
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being two dimensionless parameters ^. Parameter a varies from to 1 and 
controls the strength of self-interaction and shape of the potential. For 
a = the second minimum disappears, whereas in the limit a ^ 1 two 
minima become degenerate (see Fig. P). Parameter (5 controls size of the 
loop corrections. In order semiclassical approximation to be valid /3 should 
not be too small (see Section El for details). 

The bounce is non-trivial, 0(4)— symmetrical stationary point of 8^-, 
Eq. ()3.3p . obeying the Euler - Lagrange equation 

-r^H $ + -$^ <l>^ = , (3.6) 

dB? RdR 2 2 ' ^ ^ 

and boundary conditions 



dR 



\R=0 



0, = . (3.7) 



Here R = ((X°)^ + iXp)-*^/^. The equation ()3.(ij) at least for not very big a 
can be easily solved numerically, e.g., by the shooting method. We display 
some profiles $(-R) in Fig. |21for various values of the parameter a . 

The classical action Sci{'^) as a function of a is plotted in Fig. El (left). 
For small a classical action S goes to a constant and Sci{(y = 0) = 90.857. 
In the limit a — >■ 1 the thin-wall case is realized (see Appendix A) and 
the classical action diverges as {1 — a)~^. The ratio of the classical action 
computed numerically to the analytic thin-wall expression 



is displayed in Fig. El (right). It tends to unity for a — 1, as it should. Note, 
that the radius of the bounce increases rapidly in this limit and numerical 
calculations become delicate. So, in the present article we restrict ourselves 
to the interval a G [0, 0.95]. 
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4 Calculation of the Fluctuation Determinant 

In this section we discuss a method of computing the ratio of functional 
determinants (|2.3|) which is based on earlier papers [71 El 1101 • 

''^We use c = 1 units throughout this paper. 
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The explicit form of the operator in the nominator (j2.3p is 



M 



-A4 + + V{r) . 



(4.1) 



Here A4 is the 4-dimensional Laplace operator, and we have introduced the 
potential V as 



V{r) = U"{ip) - 



-6r](p{r) + -Av9^(r) 



m 



-3$(i?) + -a<!>\R) 
2 



= mV(i?) 



(4.2) 



The "free" operator corresponding to the metastable phase where = 

and where = U"{(f = 0) takes the same form as ()4.H) . but with V{r) = 0. 

Due to the 0(4) spherical symmetry of the bounce the operators A4 
and A^*^°) can be separated with respect to 0(4) angular momentum. We 
introduce the partial wave operators 



3d 1(1 + 2) 2 2 T./ N 
+ ^ ^ ^ + + + V{r) 

r dr 



(4.3) 



with an additional variable u that will be used later on. In terms of these 
operators we can write 



n' 

Ln 



In 



UJ 



ln{0) 



n 

/=0 



det^M;(0) 
det m{°^(0) 



(4.4) 



where di is the degeneracy of the 0(4) angular momentum, di = {I + 1)^. 
Prime denotes that for / = 1 we have to remove the four translational zero 
modes. 

The ratio of determinants of the radial operators 



det Mi{u) 
det m\''\0) 



n 



(4.5) 



can be computed using the theorem on functional determinants as described 
in the next section. Note that uf^ always denotes the eigenvalues of M/(0), 
or more generally the eigenvalues of Ai, the analogous definition holds for 
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4.1 Determinants of the Radial Operators 

In order to find Ji{iy) ()4.5|) we make use of a known theorem whose 
statement is 

= ,„„ . (4.6) 

detMS°V) 

Here ifjiiv^r) and il)f'\v,r) are solutions to equations 

MKz/)^.,, = 0. m(°)(z/)^S) = 0, (4.7) 

and have same regular behavior at r = 0. More exactly, the boundary 
conditions at r = must be chosen in such a way that the right-hand side of 
Eq. (j4.6|) tends to 1 at z/ — >• cxo. 

It is convenient to factorize the radial mode functions into the solution 
ipf'\v,r) for V{r) = and a factor 1 + hi{v,r) which takes into account 
the modification introduced by the potential. If V{r) is of finite range the 
functions 'il)f\v,r) and %l)i{v,r) have the same behavior near r = and as 
r — *■ cxD. Near r = they behave as and as r — cxd they behave as 
exp(— Kr) where k = ^/ u"^ + m?. Furthermore the requirement of analogous 
behavior near r = introduces the initial conditions /i(0) = /i'(0) = 0. The 
function h{r) then simply starts from zero at r = and goes smoothly to a 
finite constant value hi{ij, oo) as r ^ oo. The solutions il)f'\v^r) are given 
in terms of modified Bessel functions as 

^|^(z/,r) = , (4.8) 

and we have 

^,[j,,r) = [l + h{v,r)]^-^±^ . (4.9) 

r 

Then by the theorem ()4.6|) . the ratio of determinants ()4.5|) can be expressed 

as 

Ji{v) = l + hi{v,^) . (4.10) 
In terms of the h function the first equation ()4.7|) reads 

|-^+ + -fl/.,(z^,r) = y(r)[l + /iK^,r)] , (4.11) 

where I'ij^^{K,r) = dIiJ^i{Kr) / d{nr) . 
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In what follows it would be convenient to consider the perturbation ex- 
pansion 

oo 

h{v,r) = Y,hf\v,r) (4.12) 

k=l 

in powers of the potential V{r). This assumes an analogous expansion for the 
ratios Jiiv) in the sense that Ji^\v) = hf^ (z/, oo). The /c-order contribution 
h^i''^ obeys an equation 



2k^±^ — f + - 



d 
dr 



hf\v,r) = V{r)h[''-^\u,r) , (4.13) 



where we defined hf'^ = 1. Since Eq. ()4.13|) is linear differential equation 

(k) 

it holds also for linear combinations of h] . It is convenient to introduce 

notation h^j^^^ = Yl'^=k^Y^ • ■'■^ ^^^^ notation hi = h^^\ A Green function that 
gives the solution to equation ()4.13|) in the form 



hf\r) = -j dffGi{r,f)V{f)hf~^\f) (4.14) 
Jo 

with the right boundary condition at r = reads 

Gi{r,f) = h±l!^ [Ii+,{Kr^)Ki+,{Kr^) - /,+i(/tr)Ki+i(/s:f)] , (4.15) 

where r< = min{r, f}, r> = max{r, f}. 

The first term on the right-hand side of Eq. ()4.15|l does not contribute to 
h[''\oo). The Green function (I4.15p gives rise to connected graphs as well as 
disconnected ones. The latter are cancelled in ln{l + hi{oo)) whose expansion 
in /c-order connected graphs jf^loni^) reads 

In Ji{v) = ln(l + hi{u, oo)) = ^ ^-^4lni^) ■ (4.16) 

k=l 

This formula is analogous to the expansion of the full functional determi- 
nant in terms of Feynman diagrams 

lnP = ^^-^ A^''\ (4.17) 

k=l 

9 



where A^^^ is the one-loop Feynman graph of order k in the external potential 
V{r). 



Indeed, it is obvious from Eq. (j4.14j ) that h} and, therefore, J; jon ^'^^ 
of the order V'^. Since the expansion of \nV in powers of V is unique, we 
conclude that 



1=0 



2 

/ con ■ 



(4.18) 



4.2 Calculation of V^^^ 

Making use of a uniform asymptotic expansion of the modified Bessel func- 
tions in (|4.15p one can check that that j/^jon ~ 1//^^^"*^ as / ^ oo. That 
results in the expected quadratic and logarithmic ultraviolet divergences in 



InV due to the contribution of j/con Jrcon- C^ur strategy is to compute 
analytically the first two terms in the sum Eq. ()4.17|) and to add numerically 
computed hiV^^\ which is the sum without first and second order diagrams 
A^^"* and A^'^\ It reads explicitly 



(2) 



InP^^) =^(/ + l)2(lnJ,(zy)) 



(3) 



1=0 



where 



aP'(oo) - \ (ftf, 



OO 



(4.19) 



. (4.20) 



The terms in square brackets here correspond to the fish diagram Since 
all contributions to InP*^'^^ are ultraviolet finite, we need no regularization in 
computing them. The divergent contributions of the first and second order 
in V will be considered in Sec. 14.31 

In order to avoid numerical subtraction that might be delicate we re- write 
the term ()4.20j) to be summed up on the right-hand side ()4.19|) in the form 



(InJKz/)) 



(3) 



ln(l + /li(oo)) - hi{^) + -hi{(X)f 



+ hf\^)-\hf\oo)(h,{^) + hf\^)) . (4.21) 



Each of the three terms on the r.h.s. is now manifestly of order V^. The 
subtraction done in the square bracket is exact enough when the logarithm 
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is calculated with double precision. We have determined hi{r) as solutions of 

Eq.dmH) and h\^\r), /if^(r) and (r) as those of Eq. ()4.13p using Runge- 
Kutta-Nystrom integration method ^|. Of course we cannot integrate the 
differential equations until r = oo. In fact we have integrated it up to the 
maximal value for which we know the profile 0(r), and therefore V{r). This 
value is such, that the classical field has well reached its vacuum expectation 
value, and therefore V{r) has become zero. This is the condition under which 
we can impose the asymptotic boundary condition for the classical profile. 
For such values the functions hf'\r) have already become constant; indeed 
for V{r) = they have the exact form a + bKi+i^nr) / Ii+i{Kr) and the second 
part decreases exponentially for r » 1 / n. In praxi we used values of R up 

to i?max = "^^max ^ 20 - 30. 

We have neglected till now the existence of the negative mode uJq<0 
for I = and four zero modes uf = with / = 1. The former results in 
negative value of Jo{i') = 1 + hQ{i',oo) at z/ = 0. According to Eq. (j2.2p one 
has to replace Uq by \uJq\. That implies taking the absolute value of Jo(0) in 
Eq. ()4.19p : indeed Jo(0) is found to be negative. 

The translational zero modes manifest themselves by the vanishing of 
cu^Q = 0, the lowest radial excitation in the / = 1 channel with degeneracy 
= 4, and thereby by the vanishing of Ji{i^) at z/ = 0, see Eq. ()4.5j) . This 
represents a good check for both the classical solution and for the integration 
of the partial waves. The factor z/^ has to be removed according to the 
definition of det'. So in the / = 1 contribution we have to replace Ji{0) by 

Ji(z/) dJiiy) d 

Notice that replacement Eq. (j4.22|) introduces a dimension into the functional 
determinant. Thereby the units used for u become the units of the transition 
rate. Here we have used the scale m throughout, see Eqs. ()5.1|) and ()5.2p . 

Our next step is performing summation over / in Eq. ()4.19p . For small 
bounces (a < 0.8) we have found good agreement with the expected behavior, 
namely 

(lnJ/W)^oc^^-i^. (4.23) 

So, the summation has been done by cutting the sum at some value Imax and 
adding the rest sum from Imax + 1 to oo of terms fitted with 

(/ + l)5 + (/ + l)6 + (/ + l)r • ^^-^^^ 
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The summation was stopped when increasing of l^ax by unity did not change 
the result within some given accuracy S. The required accuracy was decreased 
for higher a. The problem is that the convergence becomes worse as we get 
closer to a = 1. This is related to the fact that the asymptotic behavior 
fl4.23|) sets in only when / ^ mrcs, where is the characteristic size of 
the bounce. It is of order 1/m at small values of a and can be estimated 
as 1/(1 — a)m near the thin-wall limit, a —>■ 1. As the maximal value of 
the angular momentum that we have used is / = 25, our computations cease 
to be reliable beyond a ~ 0.95. The value of 6 was about 10~^ for small 
bounces, and of order of 10~^ for a > 0.85. As we will see below, for larger 
values of a the effective action is well approximated by the leading terms of 
a gradient expansion. 



4.3 Perturbative contribution and renormalization 



We have _described in the previous subsection the computation of the finite 
part InT)^^^ which is the sum of all one- loop diagrams of the third order and 
higher, 

lnP(3) = ^Lii_AW . (4.25) 

k=3 



k 



We now have to discuss the leading divergent contributions A^^^ and A^'^\ 
These are computed as ordinary Feynman graphs. Using dimensional regu- 
larization we have 



(27r)4-^ k^ + 



where we have introduced the Fourier transform of the potential 



V{k) 



d'^xV{x)e 



-ikx 



We obtain 

^(1) 



167r2 



2 , . , 

7^; + m47r + m ■ 



+ 1 



d^xV{x) 



(4.26) 



(4.27) 



(4.28) 



where /x is the usual dimensional regularization parameter. We choose it to be 
equal to m. Then using the MS scheme we just retain the last contribution 
in the bracket (see e.g. p. 377). Thus, the finite part of A^^"* is 

1 



.(1) 



WdR V{R) 



(4.29) 
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The second order terms takes the form 



We obtain 



(2^ ^^^^ J (2^ 



(A;2 + m2)[(A; + g)2 + m2] ' 



(4.30) 



167r2 



7b + ln47r + ln-^ 



d''x{V{x)f 



(4.31) 



+ 



-/■ 

1287r4 J 



q^dq 



Viq) 



v/g^ + 4m^ , a/^^ + 4m2 + g 

2 — m = 

? V + 4m^ — g 



Again the MS" scheme corresponds to omitting the first term on the right 
hand side and for the finite part of A^"^^ we find 



(2) 



12871^ Jo 



POD 




2 


/ Q'^dQ 


v{Q) 




/o 







2 - V<3' + 4 In 



, (4.32) 



with Q = q/m being the dimensionless momenta. For the numerical eval- 
uation of A^'^^ we have to compute the Fourier transform of the external 
potential which is known numerically, the remaining computation is straight- 
forward. 



5 Numerical results 

To summarize we represented the false vacuum decay rate per unit time per 
unit volume as 



4 I Sd[^] 

7 = ( I e 



■Scl['P]-Sll_ij'p] 



where 



with perturbative 

qeff _ 1/4(1) _ 1 4(2)x 
•-^l-loop.p oV /m finJ 



and non perturbative 



2 ^ k 

k=3 



(5.1) 
(5.2) 
(5.3) 

(5.4) 
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contributions. 

It is useful to introduce the quantity G, 

G(a,/3) = ^,^{Lp[y.]/^,,H, (5.5) 

which indicates how big quantum corrections are. Since classical action, 
Eq. fl3.2p . depends linearly on parameter /3 we have G{a,P) = G{a,l)/(3. 
Numerical calculation shows that that G{a, 1) varies from 0.0367 to 0.0448 
as we vary a from to 0.95, with shallow minimum Gmin ~ 0.033 at a about 
0.6 (see Fig. Hj). Fig. El suggests that G(l, 1) ~ 0.05, which means that for 
sufficiently big values of (3, namely j3 > 0.1, the quantum corrections to the 
classical action are small (less then 50%) for all values of a. 

The corrections to the transition rate are given directly by a factor exp{—Sli^^^^p), 
so even if the classical transition rate is sizeable, as it happens for small (3, 
the quantum corrections suppress the decay of the false vacuum by factors 
exp(— 3.3) at a = and exp(— 291) at a = .9. 

Note that the main contribution to the effective action for all a is coming 
from the ^j^^, (comp. Tab. 1, Tab. 2). For small a perturbative contribution 
is almost 100% of total one loop effective action, (comp. Fig. Ej). 

In the limit a — > 1 the leading terms of the gradient expansion (Appendix 
B) gives dominant contribution to the one loop effective action. Already for 
a = 0.8 the sum of leading gradient terms 

qeff qsff I q^ff R\ 

*-'grad,0+2 ~ '-'grad.O + '-'grad,2 I'^-^J 

approximates the 1— loop effective action S'^f^^^^p within 20%. So the gradient 
expansion reproduces well the behavior of the one-loop effective action when 
a —>■ 1, see Fig. As the numerical procedure described in the main part of 
this paper becomes precarious for a > 0.9 this expansion complements the 
computation of the transition rate in this region. 

As it is well known there is exactly one negative mode in the spectrum 
of fluctuations about the bounce. Its energy is plotted vs a in Fig. IHl 

In the present paper we used dimensional regularization and have chosen 
the parameter /z^, which can be understood as parameterizing a sequence of 
possible renormalization conditions, to be equal m^. Choosing //^ differently 
would result in the following corrections to A^^^^ and A^^^ 
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where the following integral 



a(2) = / R^dR{V{R)f (5.8) 
8 Jo 

evaluated at the bounce solution. Numerical values for ^j^^, A^^^^ and a*^^) 
for different values of a are collected in Tab. 2. With the present choice of 
/i^ the perturbative terms represent the most important contributions to the 
effective action (see above), this means at the same time that a modifica- 
tion of the regularization and renormalization procedures can result in large 
changes in the one-loop effective action. 



6 Discussion and Conclusion 

In the present paper we applied previously developed technique for evalua- 
tions of functional determinants and calculated quantum corrections to the 
tunnelling transitions in 3+1 dimensional model of one self-interacting scalar 
field. 

In the present toy model decay rate is vanishingly small. The sign of 
quantum corrections is such that it decreases false vacuum decay rate. The 
corrections can be thought as originating from the particles creation during 
the phase transitions. The created particles take energy from the tunnelling 
field and therefore decrease tunnelling probability. Analytical estimations 
show that particle creation is typically weak in the thin-wall approximation 
|25j . In the present paper it was found that the quantum corrections are even 
smaller away from the thin- wall case (compare Fig. Ej), which assumes that 
particle creation for j3 > 0.1 is weak for all values of the coupling constant 
a. On the other hand for (3 < 0.1 the quantum corrections dominate, which 
means that in this regime one should look for a bounce solution taking into 
account the full effective action in the one-loop approximation fl^ |2Z] . 

Corrections to the false vacuum decay in a similar model in (3+1) dimen- 
sional theory in the thin wall approximation with the heat kernel expansion 
technique were calculated in [3^, but it is not straightforward to compare 
our results since we use a different renormalization scheme and and a dif- 
ferent parametrization of the potential. Powerful techniques for analytic 
calculations of the pre-factor using different approximations were developed 
in j^SlESEni, but we cannot compare our results directly, since these calcu- 
lations are within 3d theory. 
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The technique described here can be apphed to tunnelhng transitions in 
more reahstic theories in 4 dimensions. 
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Appendix A. The thin-wall approximation 

In the limit a 1 so called thin-wall case is realized. This is when energy 
density difference e between two vacuums 

6 = t/($_)-t/($+), (A.l) 

is small (compare to the hight of the barrier). In this case potential Eq. (|3.4j) 
can be represented as 

[/(<!>) = f/o($) + 0(e), (A.2) 
where symmetric part, Uq, in our case is 

[/o($) = ^$2(2-<|.)2, (A.3) 

o 

and 

e = 2(1 - a) , (A.4) 

In the thin- wall approximation the radius R of the bounce and the Eu- 
clidean action Sd are given analytically [H Ej as 

= ^, SIT = ^ , (A.5) 

where 

poo 

S, = 2 dR Uoi^k) , (A.6) 
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is the action of the one-dimensional kink solution corresponding to degen- 
erate potential Uq with the equal minima. For our choice of the potential, 
Eq. ()A.3|) . the kink solutions is 



One finds that Si = 2/3 and correspondingly 



1 

1 — a 3(1 — ay 



Appendix B. The leading terms of the gradient 
expansion 

We want to derive an approximation to the effective action of a scalar field 
on the background of a bounce solution. The strategy is to expand first the 
effective action with respect to external vertices, and to expand in a second 
step the resulting Feynman amplitudes with respect to the external momenta. 
This approach is fairly standard, and has been used, e.g., in Ref. [21]. We 
note that we will retain all powers in the external vertices; such a summation 
was found to yield a very good approximation for the sphaleron determinant 
[m 112], see Fig. 1 in the second entry of Ref. ^21- We have to compute 
the trace log or log det of a generalized Euclidean Klein-Gordon operator 
A4 + U"{(j)) where A4 is the four-dimensional Laplace operator. Formally 

[\nV] = In 

We introduce a potential V{x) via 

= + V{x) ; f/"(0) = . (B.2) 

For the bounce the potential depends only on r = |x| but we will not use this 
now. The logarithm can be expanded with respect to the potential V{x). 



-A4 + f/"(0) 



(B.l) 
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We write 



In 



-A4 + m^ + V{x)' 
-A4 + m2 



= ln[(-A4 + m2)-i(A4 + m2 + y(x))] 
= In [1 + (-A4 + m2)-V(x)] 
" -1)^+1 



(B.3) 



E 

Ar=i 



AT 



[(-A4 + m2)-V(x)] 



AT 



and the effective action is given by 

" .l)iV+l 



5 



eff 



Ar=l 



2A^ 



tr [(-A4 + m2)-V(x)] 



AT 



We introduce the Fourier transform 



V{q) 



(B.4) 



(B.5) 



The individual terms in the expansion of the effective action have the form 
of Feynman diagrams with external sources V{qj) with j — 1 . . .k. The 
momentum that has flown into the line I is 

I 

Qi = Y.^j. (B.6) 

of course the total momentum must be zero,i.e., Qj^ — 0. With these nota- 
tions we can write the A'^th term in the effective action, omitting the factor 
(-l)^+72Ar as 



A 



N 



N 

n 



n N 



(27r)^ 



n 



(27^)^5(g^^) ■ 



(B.7) 

The four-momentum delta function arises from taking the trace. We obtain 
a gradient expansion by expanding the denominators (p + Qi)"^ + w? with 
respect to the momenta Qi. The leading term is of course 

hAT r 



^Af,0 




d Qj ~ 



d^x [V{x)] 



N 



(B.8) 
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The zero-gradient contribution to the effective action is obtained by resuming 
this series; one finds 



s: 



■eflf 
;rad,0 



Of course this integral has to be regularized, e.g., via dimensional regulariza- 
tion. The divergences come form the terms with = 1 and N = 2, which 
are standard divergent one loop integrals. 
We find 



r(f) 



^ fdpp^ 



p2 + c/"(0) 



r (f ) (47r)^/2 
2 



dpp 



D+l 



p^ + U"{(p) 

1 



p=0 

1 



y + u"{(p) p2 + c/"(o)_ 

The first term in the parenthesis vanishes for < D < 2 and is defined to 
vanish in general by analytic continuation. The second term can be rewritten 
as 



{[u"{4>)f/^-[u"{o)f/^} 



dx 



X 



D+l 



X^ + 1 



DF (f ) (47r)^/2 

. r(D/2 + i)r(-D/2) rr^.. _ ^u"(o)rn 



(47r)^/2 
Now set = 4 — e and use 



1 



(-2 + e/2)(-l + e/2) V2 

1 r2 



- <^ - - 7e + 



to obtain 



7fi + In 47r + - 



327r2 , 

X j(m2 + V(r))2 



e, m2 + T/(r), 

1 In „ 

2 ; 



— m 



l--ln^ 
2 //2 
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Using MS subtraction we get 



■1 r3 



- [2mV(r) + V^2(r)] 



327r2 1^2 

Integrating over 4d Euclidean space we finally obtain 



(B.IO) 



qeff 



— / R^dR 
32 70 



l + y(i.)flnl±^ 



{2V{R) + V^{R))+\njl^ 



(B.ll) 



with jl = fi/m. 

Let us now consider the one- and two-gradient contributions. We expand 
the denominators up to second order in the gradients, i.e., in the momenta 
Qj. We obtain 



N 



- n 



/=1 



1 



1 



1 



N 



{p + QiY + 

N 



(p^ -|- m?)^ {p^ + m?) 



N+l 



N-1 N 



^ j=i ^-f^ ^ j=i k=j+i 



N 



(p2 + m2)^+2 ^ 
^ ' .7=1 



(B.12) 



Under 0(4)— symmetric integration -ipf^Pu — p'^S^i,, and p^ ^ 0. So the one- 
gradient term vanishes and the complete two-gradient contribution becomes 



n 



iV,2 



(p2 + m2)^+2 
1 

(p2 + m2)iV+2 



3 k>j j 



k>j j 



(B.13) 



We now have to rewrite this in terms of the momenta qj that represent the 
gradients on the functions V{qj). After having used the fact that 112 appears 
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under the integral over d'^p wc will now use the fact that it appears under the 
product of integrals J d^qjV{qj) which implies permutation symmetry in the 
indices j. So if we expand the products Qj ■ Qk and we will encounter just 
two kinds of terms: products qi ■ qm with / ^ m and squares q^, which may be 
replaced by q\ • q^ and by qf, respectively. We have to do some combinatorics 
in order to find 

{N-1)N{N+1) , N(N + 1) 2 
Z^Qj - V ^?i-52 + ^^ -ql (B.14) 



3 - - ■ 2 

{N -1)N{N +1){3N -2) 
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k>j 



-qi ■ q2 



(N -1)N(N +1) . 
+- '-^ -qf ■ (B.15) 

Now we may use momentum conservation to rewrite 

ql = -gi • (^2 + • • • + qn) ^ -{N - l)gi • q2 (B.16) 



so that 



^^2 {N-1)N{N + 1) 

l^Qj ^ ^51 ■^2 (B.17) 



k>j 



6 

^ (N -2)(N -1)N(N + 1) 
J^QrQk ^ ^ -qi-q2 (B.18) 



and 



The momentum integrals are 



(27r)4 (p2 + ^2)7v+2 i67r2 {N -1)N{N+1) ^ ' 

(27r)4 (p2 + m2)^+2 167r2 7V(Ar + 1) ^ " ^ 
and, therefore, 

--4-n^2 = gi • q2 . B.22 
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/ 



The momenta are converted into gradients; so we finally obtain as the ex- 
pansion terms of the two-gradient part of the effective action 



Amo — 



167r2 



1 d^'x 


'v{xy 







N-2 



N 



12m2 



(B.23) 



The term ^41^2 is zero. The sum over all terms yields 



1 



grad,2 33^2 



d X 



1 



1 



m2 + V{x) 12 



{VV{x)f 



(B.24) 



or finally in dimensionless variables 



ceff ^ _}_ 
'^grad,2 j^gg 



WdR 



1 



1 + V{R) 



{V'{R)f 



(B.25) 



An alternative derivation starts with a technical step that frees us from 
the denominator We take the derivative of the effective action with 

respect to m^, a step that we can revert later on. We then obtain, using the 
cyclic property of the trace, 



Q = 



dS' 



eS 



dm? 



00 / -, \ AT 

= ^^^2^tr{[(-A4 + m2)-V(x)]^(-A4 + m2)-i} (B.26) 



7V=0 

00 



Ar=o 



We note that we have included the = term, which can be removed later 
on if necessary. So we have arrived at the trace of the exact Green function 
in the external field. The terms i^jv have the form 



B 



N 



'-^'"/(Sn[/ 



X 1/(^2 



Viqs) 



p2 _|_ 



V{qi 



{p + Q2Y + m? 
[P + Qn) + rn^ 



(B.27) 
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Assume we have expanded the fraction + Qk)"^ + m^] to first order in 
'^V ' Qk + Q\) yielding a factor 

1 r . „ 1 



p2 _j_ ^2 



[-2p ■ Qk - Ql] 



(B.28) 



at the k place in the product of propagators and vertices, in other words 
we have obtained an insertion of —2p ■ Qk — Q\- Consider the part of the 
product to the right of this insertion. We rewrite it as 



n 



N 

j=k+l 



X n,-=fe+i 



d^XjV{x 



[-2p ■ Qk - Ql] 



p^ + 



p2 _|_ ^2 

{27cyS{Qk + Qk+i + ■■■ + qNlB.29) 



We furthermore rewrite the delta function as 



Inserting this in (|B.29j) we can carry out the integrations over the qj and the 



Xj to obtain 



1 ^ 

/ d'xe'^"-^ [-2p ■ Q, - Ql] ^— ^ J] 



j=k+i 



V{x] 



p2 _|_ ^2 



(B.31) 



idn on the 



Now the Qfc^^ in 2p ■ Q^ + Ql can be written as —id/dx^ 
exponential. Integrating by parts they can be written as id^j, acting on the 
product to their right. So the whole string to the right of the insertion can 
be written as 

N 



1 



p2 _|_ ^2 



n 

j=k+i 



V{x] 



p2 _|_ ^2 



(B.32) 



We now consider the sum over N; we split = k+l and (—1)^ = (— 1)'^(— 1)'. 
The sum over / is independent of k and runs from to oo and, putting in 
the factor (—1)' we obtain 



ld^xe-^-^[-2.p.8,8^] 



-Vix] 



1 



p2 _|_ j^2 



p2 _|_ ^2 



[B.33) 
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Note that the sum starts with I — 0, which corresponds to the case k = N; 
in this case the product over j reduces to 1. Now we do the analogous 
operations on the part to the left of the insertion, using in the exponent 
Qfc = Q'l + • • • + Q'fc; we now can carry out the summation over k and we find 
finally for the case that we have taken into account the first order expansion 
of one of the denominators {p + Qk^ + rn^ 



I 



p2 _|_ ^2 _|_ Y(^x) '- + + V{x) 



Obviously part —i2p ■ d vanishes upon symmetric integration over p. It also 
can be written as a boundary term for the x integration. If we want to obtain 
the second order gradient term we have to take into account the Ql term of 
the first order expansion, i.e. 

^4 \ q2 \ (-r, oe.\ 

p2 + ^2 + v{x) p^ + m'' + V{x) ' ^ ' 

the terms —2ip ■ d arising if two denominators are expanded to first order, 
yielding 

4 1 1 



p2 _|_ j^2 _|_ V(x) p^ + nn? + V{x) p^ + nn? + V{x) 

(B.36) 

Here is included the term arising from expanding one propagator to second 
order. Indeed this yields 

[-2p • Qk - Ql] [-2p • Qk - Ql] , (B.37) 

a term that is needed for obtaining the complete propagator l/{p^ + + 
V{x)) between the two insertions. We now have the two-gradient term 

2 J (27r)4y \p2 + ^2 + ^(^) p2 + ^2 + ^(^) ^ • 

+ , , ^ , (-2ip • d)^^-^^—--(-2vp -d)- ^ 



p2 _|_ ^2 V{x) p^ + m? + V{x) p^ + + V{x) 

The first term can be written, after one integration by parts as 



1 r _d^ r 

2 7 (27r)4y 



(^^a;-^^ :T^^-^7r^[dV(x)? . (B.39) 
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In the second term we remark that the derivatives in the first insertion act 
on the complete part to the right of it. Therefore an integration by parts lets 
it act onto the part to the left of it. Using symmetric integration over p the 
second part yields 



d X 



p^ 



[p^ + m? + V{x))^ 



[dV{x)\ 



(B.40) 



Now we integrate with respect to m? to obtain the two-gradient contribution 
to the one-loop effective action 



ceff 
^grad,2 



327r2 



d'^p 
d'^x 



dS 



1 



1 



p 



3(p2 + m2 + V^(x))3 

1 2 



4 (p2 + w? + V{x)Y 



m2 + V{x) 12 



[dV{x)X 



[dV{x)f 
(B.41) 



which coincides with the previous result Eq. ()B.24|) . 

The terms of the gradient expansion can be evaluated in a straightforward 
way. We note, however, that the term m? + V{x) vanishes, depending on 
value of a, at one or two points, and that therefore the expressions are ill- 
defined a priori. This is a reflection of the fact that the effective action 
has an imaginary part, due to the negative mode. An expansion of the 
effective action has to reflect this feature. With an m? — it prescription this 
becomes apparent. When computing these terms we have used the principal 
value prescription for 2 taken the absolute value in the logarithm 
appearing in S^^l^ ^. 
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Figure 1: Potential in dimensionless form Eq. ()H.4j) . The curves are 

labelled with the value of a. 
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Figure 2: Bounce profiles for different a. 




Figure 3: Classical action Sd versus a (left) and the ratio Sd/S^ for a > 0.5 
(right). 
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Figure 5: Our results for the effective action Sl[_{ (squares) togetlier with 



the perturbative part Slil^^^ ^ (dotted hue) and the leading parts of the gra- 

ed li] 

are multiplied by the factor (1 — a)^. 



dient expansion S^^f^ Q^2 (dashed line, a = 0.45 — 0.95). All shown quantities 
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Figure 6: The negative mode energy as a function of a 



33 



a 


qeff 
1— loop,p 


qeff 
1— loop,n.p. 


qeff 




0.00 


3.253E+00 


8.216E-02 
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1.153E+02 


0.20 
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0.25 
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4.865E+00 
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0.30 
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-3.201E-01 


5.351EH-00 


1.552E+02 


0.35 


6.499E+00 


-5.539E-01 
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1.744E+02 


0.40 


7.571E+00 


-8.836E-01 


6.687E+00 


1.983E+02 


0.45 


8.984E+00 


-1.348E+00 


7.637E+00 


2.286E+02 


0.50 


1.089E+01 


-2.006E+00 


8.889E+00 


2.681E+02 


0.55 


1.356E+01 


-2.958E+00 


1.060E+01 


3.211E+02 


0.60 


1.741E+01 


-4.371E+00 


1.303E+01 


3.951E+02 


0.65 


2.326E+01 


-6.560E+00 


1.670E+01 


5.033E+02 


0.70 


3.277E+01 


-1.015E+01 


2.261E+01 


6.720E+02 


0.75 


4.966E+01 


-1.659E+01 


3.306E+01 


9.589E+02 


0.80 


8.382E+01 


-2.969E+01 


5.413E+01 


1.512E+03 


0.83 


1.240E+02 


-4.512E+01 


7.887E+01 


2.136E+03 


0.85 


1.686E+02 


-6.233E+01 


1.062E+02 


2.809E+03 


0.87 


2.409b+02 


-9.038E+01 


1.50ob+02 


3.874E+03 


0.88 


2.950E+02 


-1.114E+02 


1.836E+02 
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4.711E+02 
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0.91 
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-2.390E+02 
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0.92 


8.455E+02 


-3.284E+02 


5.171E+02 


1.227E+04 


0.93 
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-4.724E+02 


7.347E+02 


1.711E+04 
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1.829E+03 


-7.209E+02 
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-1.188E+03 


1.820E+03 


4.061E+04 



Tab. 1. Numerical results for classical action and one loop effective action. 
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.(1) 


J (2) 




0.00 


5.178E+00 


-2.654E+00 


1.036E+01 
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5.706E+00 
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1.085E+01 


0.10 


6.325E+00 


-2.358E+00 


1.144E+01 


0.15 


7.060E+00 


-2.235E+00 


1.215E+01 


0.20 


7.942E+00 


-2.133E+00 


1.300E+01 


0.25 


9.013E+00 


-2.059E+00 


1.405E+01 


0.30 


1.033E+01 


-2.020E+00 


1.536E+01 


0.35 


1.199E+01 


-2.024E+00 


1.702E+01 


0.40 


1.410E+01 


-2.085E+00 


1.916E+01 


0.45 


1.686E+01 


-2.219E+00 


2.199E+01 


0.50 


2.056E+01 


-2.453E+00 


2.585E+01 


0.55 


2.570E+01 


-2.825E+00 


3.127E+01 


0.60 


3.311E+01 


-3.399E+00 


3.921E+01 


0.65 


4.438E+01 


-4.286E+00 


5.147E+01 


0.70 


6.269E+01 


-5.701E+00 


7.175E+01 


0.75 


9.526E+01 


-8.109E+00 


1.085E+02 


0.80 


1.613E+02 


-1.269E+01 


1.848E+02 


0.83 


2.391E+02 


-1.778E+01 


2.761E+02 


0.85 


3.256E+02 


-2.321E+01 


3.790E+02 


0.87 


4.oo0b+02 


-3.170b+01 


5.479b+02 


0.88 


5.711E+02 


-3.788E+01 


6.753E+02 


0.89 


7.137E+02 


-4.609E+01 


8.493E+02 


0.90 


9.134E+02 


-5.735E+01 


1.094E+03 


0.91 


1.203E+03 


-7.333E+01 


1.453E+03 


0.92 


1.643E+03 


-9.699E+01 


1.999E+03 


0.93 


2.347E+03 


-1.341E+02 


2.883E+03 


0.94 


3.561E+03 


-1.963E+02 


4.417E+03 


0.95 


5.861E+03 


-3.118E+02 


7.359E+03 



Tab. 2. Numerical results for the first and second order contribution 

coefficients. 
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